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ABSTRACT 

Although at least one quarter of early-type barred galaxies host secondary stellar bars embedded in their 
large-scale primary counterparts, the dynamics of such double barred galaxies are still not well understood. 
Recently we reported success at simulating such systems in a repeatable way in collisionless systems. In order 
to further our understanding of double-barred galaxies, here we characterize the density and kinematics of the 
V-body simulations of these galaxies. This will facilitate comparison with observations and lead to a better 
understanding of the observed double-barred galaxies. We find the shape and size of our simulated secondary 
bars are quite reasonable compared to the observed ones. We demonstrate that an authentic decoupled sec- 
ondary bar may produce only a weak twist of the kinematic minor axis in the stellar velocity field, due to the 
relatively large random motion of stars in the central region. We also find that the edge-on nuclear bars are 
probably not related to boxy peanut-shaped bulges which are most likely to be edge-on primary large-scale 
bars. Finally we demonstrate that the non-rigid rotation of the secondary bar causes its pattern speed not to 
be derived with great accuracy using the Tremaine- Weinberg method. We also compare with observations of 
NGC 2950, a prototypical double-barred early-type galaxy, which suggest that the nuclear bar may be rotating 
in the opposite sense as the primary. 

Subject headings: stellar dynamics — galaxies: evolution — galaxies: kinematics and dynamics — galaxies: 
structure 
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1. INTRODUCTION 

Double barred (S2B) galaxies were first described over 
thirty years ago (Ide Va ucouleurs 1975). The Hubble Space 
Telescope has revealed secondary bars at the center of 
at least one quarter of early-type optically-barred galaxies 
(lErwin & Sparkd [2b02|) . Dynamically decoupled^ secondary 
bars in S2B galaxies have been hypothesized to be a mecha- 
nism for driving gas past the inner Lindblad resonance (ILR) 
of primary bars, to feed th e supermassive black holes that 
power active galactic nuclei (IShlosman et al.lfT989|) . 

The dynamics of secondary bars are still not well 
understood. The random apparent relative orientations 
of primary and secondary bars in nearly face-on galax- 



ies points to dynamical decoupling (iButa & Crocked [19931: 
iFriedli &^ar tinet 1993). But images alone cannot re- 
veal much about how the two bars rotate through each 
other. Kinematic evidence of decoupling, using either 
gas or stars, is harde r to obtain jP etitpas & WilsonI l20Q2l: 
ISchinnerer et al.l l2QQ2t iMoiseev et al . 2004). Indirect evi- 
dence for decoupling was claimed by Ems ellem et al.l (l200l|) 
based on rotation velocity peaks inside the secondary bars. 
Conclusive direct kinematic evidence for a d ecoupled sec- 
ondary bar was obtained for NGC 2950 by ICorsini et al.l 
(I2003L her eafter CDA03) who showed, using the method of 
iTremaine & Weinberg (1984) , that the primary and secondary 

^ Harlan J. Smith Fellow 
Electronic address: shen@astro.as.utexas.edu 

^ Brooks Prize Fellow; current address: Centre For Astrophysics, Univer- 
sity of Central Lancashire, Preston, UK PRl 2HE 
Electronic address: debattis@astro.washington.edu 

^ In this context, by decoupled we mean only that Q^s^^p^ where Q^s (^p) 
is the pattern speed of the secondary (primary) bar. 



bars cannot be rotating at the same pattern speed. 

Simulations offer the best way to understand double barred 
systems. However, the decoupled nuclear bars that formed 
in early simulations did not last long. For example, the most 
long-lived nuclear bar in Friedli & Martinet (1993) lasted for 
less than two turns of the primary bar, corresponding to about 
0.4 Gyr, which is far too short to explain the observed abun- 
dance of nested bars. Furthermore, their models usually re- 
quire substan tial amounts of gas to form and maintain these 
nuclear bars. iHeller et aP (|2007a,b) reported that nested bars 
form in a quasi-cosmological setting, but the amplitudes of the 
bars also seem to weaken rapidly af ter most of gas has formed 
stars (iHeller et al.ll2007aL Figure 2). iPetitpas & WilsonI (l2004l) 
found that 4 out of 10 double-barred galaxies contain very 
little molecular gas in the nuclear region. These clues sug- 
gest that large amounts of molecular gas may not be neces- 
sary to maintain central nuclear bars. Rautiainen et all (|2002|) 
reported that a secondary bar forms in a collisionless A/^-body 
simulation, although their secondary bar had a "vaguely spiral 
shape." 

On the side of orbital studies, Maciejew ski & Sparkd 
(1 19971 i2000) discovered a family of loop orbits that may 
form building blocks of long-lived nuclear stellar bars (also 
'Maciejews ki & Athanassoul al [2007b . Their studies are very 
important for the understanding of double barred galaxies, 
but their models are not fully self-consistent, since nested 
bars in general cannot rotate rigidly through each other 
(iLouis & Gerhard! Il988h . So fully self-consistent V-body 
simulations are still needed to check if their main results still 
hold when the non-rigid nature of the bars is taken into ac- 
count. 

Recently iDebattista & ShenI (I2007L hereafter DS07) 
demonstrated that long-lived secondary bars can form in 
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Fig. 1. — 3-D surface density of 3 runs. In each row, from left to right, panels are for the large-scale total stellar distribution, central zoom-in of total stellar 
distribution, separate disk (dark) and bulge (red) surface density contours, respectively, (a), top row: the canonical run with a secondary bar formed (Run D); 
(b). middle row: the run with an unrotating bulge, which formed only a single bar (run S); (c). bottom row: the run without a bulge (run NB). The contours in x-y 
plane are separated by half dex, while in vertical projections they are spaced more sparsely to avoid contour over-crowding. 



purely collisionless A/^-body simulations, when a rotating 
pseudobulge is introduced in their model. The nuclear bars 
in their work are distinctly bars, and do not have a spiral 
shape. They showed that the behavior of their models 
were in g ood agreement with the loop orbit predictions of 
iMaciejew ski & Sparke (2000). 

In this report we analyze the photometrical and kinematical 
properties of high resolution models in detail. Our theoret- 
ical results here can also be compared to the observed 2-D 
kinematics of some double-barred galaxies, to achieve a bet- 



ter understanding of the dynamics of the secondary bars. 
2. MODELS 

The simulations presented in this paper are all collisionless. 
The model setup is very similar to that of DS07. As in DS07, 
the formation of the secondary bar is induced by a rotating 
pseudobulge. We focus on three simulations: run D which 
formed a long-lasting double-barred system due to an initially 
rotating bulge, run S in which only a single bar formed with an 
initially unrotating bulge component, and run NB where there 
is no bulge component initially. Our high-resolution simula- 
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tions consist of live disk and bulge components in a rigid halo 
potential. We restrict ourselves to rigid halos to allow higher 
mass resolution in the nuclear regions, to study the compli- 
cated co-evolution of the two bars without the additional evo- 
lution introduced by the halo. The rigid halos used in this 
study are all logarithmic potentials <l>(r) = jV^ ln(r^ + r^). We 
set Vh = 0.6 and rh = 15 in all runs. We employed about four 
times more particles than the runs published in DS07 to better 
analyze the photometric and kinematic properties; Run D and 
S have 4.8 x 10^ equal mass particles, with 4 x 10^ in the disk 
and the rest in the bulge. Run NB has 4 x 10^ in the disk only 
since there is no bulge. 

The initial disks in our simulations all have exponen- 
tial surface densities with scale-length R^, mass and 
Toomre-2 2. The bulge w as gen erated using the 
method of iPrendergast & Tomeil (ll97Ql) as described in 
iDebattista & SellwoodI (|200Q|) . where a distribution function 
is integrated iteratively in the global potential, until conver- 
gence. In both run D and run S the bulge has mass Mb = 0.2M^ 
and we used an isotropic King model distribution function. 
The bulge truncation radius is 0.9Rd in both run D and run 
S. The bulge set up this way is non-rotating. We introduce 
bulge rotation in run D by simply reversing the velocities 
of bulge particles with negative angular momenta, which is 
still a vali d solution of the collisionless Boltzmann equation 
(iLynde n-Bell 1962). The bulges in run D and run S are flat- 
tened by the disk potential initially, and remain so at later 
times. The initial kinematic ratio Vp/a i n run D is slig htly 
above the line for oblate isotropic rotators (iBinneyll 19781) 

We use 7?d and as the units of length and mass, respec- 
tively, and the time unit is (Rl/GM^)^^^. If we scale these 
units to the physical values Md = 2.3 x 10 ^^M© and Rd = 2.5 
kpc, then a unit of time is 12.3 Myr. We use a force resolution 
(softening) of 0.01, which scaled to the above physical units 
corresponds to 25 pc. These sim ulations were evolved wit h 
a 3-D cylindrical polar grid code (|Sellwood & Vallurilll997|) . 
This code expands the potential in a Fourier series in the cylin- 
drical polar angle ^; we truncated the expansion at m = 8. 
Forces in the radial direction are solved for by direct convo- 
lution with the Greens function while the vertical forces are 
obtained by fast Fourier transform. We used grids measuring 
A^/? X X A^^ = 58 X 64 x 375. The vertical spacing of the 
grid planes was = 0.01/?d. Time integration used a leapfrog 
integrator with a fixed time step 5^ = 0.04. 

3. PHOTOMETRY 

3.1. Shape of the secondary bar 

Figure [T] shows the surface density contours and images of 
the double-barred run D, run S, and run NB which does not 
have an initial bulge component. From Figure \T\ we see that 
for run D the secondary bar shows up in both the disk and 
bulge components. All three large scale bars appear qualita- 
tively similar to each other in both the face-on and edge-on 
views. 

When viewed side-on, the large scale bar appears to 
be boxy /peanut- shaped regardless of whether or not an 
initial live bulge is included, or if a secondary bar 
is present. The formation of a boxy/peanut- shaped 
bulge from disk has been studied extensively wit h A^- 
body simulations (e.g. [C ombes et al. 1990; Raha et al] 
1991; "Bureau & Athanassoula 2005; Debattista et al. 20051: 
Martinez- Valpuesta et al. 200 6: Debattista et al. 2006). Fig- 
ure [T] does not show any obvious influence of the sec- 
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Fig. 2. — Run D at r = 405 projected to / = 45° and ^nuc =45° with 
all particles shown. Th e model bears a passing resemblance to NGC 2950 
(Erwin & Spa rkd[2003) . 

ondary bar on the overall boxy /peanut- shaped side-on ap- 
pearance of a large-scale bar. So it is quite unlikely 
that most boxy- s haped bulges are edge - on nu clear bars, 
as speculated by iKormendy & KennicuttI (l2004l) as a pos- 
sibility of explaining boxy bulges. This is hardly sur- 
prising as the small size of the secondary bar makes its 
side-on signatures, if any, easily masked by the primary. 
A caveat may be that the boxines s in Figure [D does no t 
cover a range as wide as that in IDebattista et al.l (l2006h . 
Also note that the boxy part is smaller than the pri- 
mary bar (regardless of double-barred or single-barred) as a 
whole (see IShen & Sellwood 2004; Kormendv & Kennicutt; 
"2004"; "Martinez-Vabuesta et al., .2006: .Debattista et al...2006: 
Athanassoula & Beaton 2006). 

Figure [21 shows the projected system (at t = 405 when the 
two bars are nearly perpendicular) with an ordinary orien- 
tation: the system is inclined at / = 45° with the line of 
nodes (LON) of ipnuc = 45° relative to the secondary bar ma- 
jor axis. The surface density image and contours resemble 
many observed double-barred systems, such as NGC 2950, 
even though we did not deliberately set out to match it. 

3.2. Size relation of the two bars 

Figure [3] shows radial variations of m = 2 Fourier amplitude 
and phase for run D at ^ = 400. Figure (H shows the ellipticity 
and position angle (PA) profiles of ellipses fitted with IRAF 
for the same data as in Figured (we use log scale for radius to 
be consistent with what observers usually adopt). There are 
four popular methods for determining the semi-ni ajor axis as 
of a ba r. As summarized by O'Neill & Dubinskil (|2003) and 
lErwinI (1200 5). For convenience, we denote the primary bar as 
Bl and the secondary bar as B2. 

(1) the bar end is measured by extrapolating half-way down 
the slope on the m = 2 amplitude plot (Fig [3^). We find aBi ^ 
2.3, aBi ^ 0.4, the B2/B1 bar length ratio is about ^ 0.17. 

(2) the bar end is measured when m = 2 phase deviates from 
a constant by 10° (Fig[3]3). We find ani ^ 2.1, asi ^ 0.4, the 
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Fig. 3. — The radial variations of the m =2 Fourier ampUtude and phase of 
all particles for Run D at f = 400. 



B2/B1 bar length ratio is about 0.19. 

(3) the bar end is measur ed at the peak of the 
fitted eUipticity profiles (e-g- iMarinova & Joged 120071: 
iMenende z-Delmestre et al.ll2007|) . which is shown in Fig|4^. 
We find gbi ^ 1.7, asi ^ 0.2, the B2/B1 bar length ratio is 
about-- 0.12. 

(4) the bar end is measured when the PA of fitted ellipses 
deviates from a constant by 10°. (Figllt). We find gbi ~ 2.3, 
aB2 ^ 0.4, the B2/B1 bar length ratio is about 0.17. 

Method 1 , 2 and 4 yield consistent values of the bar lengths 
and length ratios. We found that method 3 tends to give a 
lower val ue of bar lengths than the other three methods, as 
shown in lQ'Neill & Dubinskil (l2003h . Although these meth- 
ods have some uncertainties in measuring the bar lengths, the 
length ratio of the two bars is in the range of 0.12 to 0.19 (in 
particular method 1, 2, and 4 give a consistent narrow range 
of 0. 17 to 0. 19). This result is in good agreement with the typ- 
ical observe d length ratio of local S2B systems (median ratio 
- 0.12, see 'Erwin & Sparkd 120021: lErwinI 120041: iLisker et al.l 
^006). Note that we expect that the length of the secondary 
cannot be too large, otherwise the gravitational torque from 
the primary bar will inevitably twist the secondary into align- 
ment if they rotate at different pattern speeds. 

3.3. Face-on surface density profiles 

Figure [5] shows the face-on surface density profiles along 
the major and minor axis of the primary bar in run D, S, and 
NB. Compared with the initial profile of each run, there is a 
significant increase in central density following the formation 
of the large scale bar which redistributes the disk particles. 
For run D, the central density profile along the primary bar 
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Fig. 4. — EUipticity and position angle as a function of semi-major axis of 
IRAF-fitted ellipses for Run D at f = 400. 

major axis is no longer higher than that along the minor axis, 
due to the secondary bar orienting to a different direction from 
the primary bar. This minor axis over-density is, of course, 
even more pronounced when the two bars are perpendicular. 
This can be an important signature of confirming small sec- 
ondary bars photometrically, especially when the central re- 
gion is not well resolved. Variations in M/L are unlikely to 
mask this minor/major axis difference as density profiles are 
for roughly the same radial range. On the other hand, long 
wavelength photometry is preferred to minimize the effects 
of dust. We do not find other significant differences in the 
face-on surface density profiles between run D and other runs 
without a secondary bar. 

4. KINEMATICS 

Figure [6] shows the behavior of the azimuthally averaged 
Vt, and the location of the Lindblad resonances of 

the bars at around t = 400 for run D. As shown in DS07, 
the pattern speeds of the bars, especially that of the sec- 
ondary, vary as they rotate through each other: the sec- 
ondary bar rotates slower than average when the two bars 
are perpendicular, and faster when the bars are parallel. The 
patten speed bands shown in Figure [6] reflect such varia- 
tions. Clearly the pattern speed of the secondary bar os- 
cillates much more than that of the primary. The primary 
bar extends roughly to its CR radius (~ 2.5), consistent with 
the general expectati on and is therefore considered a fast bar 
(e.g. CDA03, Debat tista & W illiams 2004). The secondary 
bar rotates faster than the primary bar. However, the sec- 
ondary bar is much shorter than its shortest Rqr- In ad- 
dition, even if the variation of the pattern speed is taken 





Fig. 5. — Face-on surface density profiles along major and minor axis of t 
For each row, from left to right: the surface density profile att = 0,t = 400, ai 
lines are surface density for all particles, dashed lines are for bulge particles c 
and minor axis of the primary bar, respectively. 

into account, the Rqr of the secondary is not very close to 
the 7?iLR of the primar y, if we use the same naive defini- 
tion of 7?iLR as in Pfen niger & NormanI (|i990)'^. This is in- 
consistent with the speculated C R-ILR coupling requirement 
for making secondary b ars (e.g. iPfenniger & NormanI 119901: 
iFriedh & Martinetlll993h . 

Since the long-lasting secondary bars which form in our 
simulations generally do not extend to their corotation radii, 
the kinematics of secondary bars differ from those of primary 
bars in at least this important detail. We here explore the kine- 
matic observables of S2Bs in our simulations in more detail. 

4.1. Line-of-sight velocity distribution 

We analyzed the line-of-sight velocity distribution 
(LOSVD) by measuring the mean velocity v and velocity 
dispersion cr. Departures from a Gaussian distribution are 
parameterized by Gau ss-Hermite moments (lGerhardlll993l: 
Ivander Marel & Franxl 119931: iBender et al.1 11994 . The 
second order term i n such an expa nsion is related to the 
dispersion. Following iGerhardI (Il993|) the third-order term /z3 
and fourth-order term /z4 are defined as 

^ A cautionary note is that the 7?ilr read naively from Figure |6] serves 
just as a visual guide, because the /?ilr determined this way is reli- 
able only for weak bars, w hich is questionable for our strong bars (e.g., 
Ivan Albada & SanderslfT982h . 



primary bar for run D (top row), run S (middle row) and run NB (bottom row), 
the close-up view of the inner region at r = 400, respectively. In all figures, solid 
y, dotted are for disk particles only. The black and red curves are along the major 
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Fig. 6. — Frequencies as a function of radius at around t = 400 for run 
D, calculated based on the azimuthally averaged gravitational attraction. The 
full-drawn line shows the curve of the circular angular frequency Q and the 
dashed curves mark O ± where k is the epicyclic frequency. The two 
shaded bands show the oscillational ranges of the bar pattern speeds (the up- 
per band is for the secondary bar and the lower one is for the primary). 

where w = (v-v)/a, n = 3 or 4, Hsiw) = [l/(967ry/^](^w^ - 
12w), and H4(w) = [l/(16S7ry/\l6w'^ -4Sw^ ^12). For a 
particle model, the integral becomes a sum and S is replaced 
by Np, the number of particles in a bin. hs measures deviations 
that are asymmetric about the mean, while h4 measures the 
lowest order symmetric deviations from Gaussian (negative 
for a "flat-top" distribution, and positive for a more peaked 
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one). 

Figure|7la-e) show the images and the LOS stellar kinemat- 
ics of run D at ^ = 405 when the two bars are almost perpen- 
dicular to each other. For comparison purpose Figure |7lf) is 
for run NB (run S is very similar to run NB, so it is not shown 
for brevity). As in Figure O we project the system to / = 45° 
with the LON of 45° relative to the secondary bar major axis. 

The most striking feature in Figure |7]is that the twist of the 
kinematic minor axis (i.e. vios = 0) in the secondary bar region 
is weak (see the mean velocity maps). The kinematic minor 
axis is almost perpendicular to the inclination axis, although 
there is a small but noticeable twisted pinch near the kine- 
matic minor axis in the nuclear region. The weak central twist 
is mainly due to the relatively large velocity dispersion, espe- 
cially in the central region (likewise at ^ = 20 when only the 
small nuclear bar exists, the stellar twist is slightly stronger 
than at ^ = 405, but still quite small compared to the expected 
twist in gaseous kinematics). On the other hand, the twist 
of the kinemati c majo r axis i s more prominent in the central 
region. Moiseev et al. (2004) found the stellar kinematic mi- 
nor axis hardly twists from the PA of the disk in their sample 
with the most reliable kinematics, leading them to question 
whether nuclear photometric isophotal twists represent bona 
fide dynamically decoupled secondary bars. We demonstrate 
that an authentic decoupled secondary bar may indeed pro- 
duce a very weak twist of the kinematic minor axis in the 
stellar velocity fiel d. So a central stellar velocity map without 
a strong twist as in lMoiseev et al.l (12004) does not necessarily 
exclude the existence of a decoupled nuclear bar. 

As a comparison. Figure |7lf) shows that the kinematical 
minor axis twist is just slightly stronger for the single-barred 
run NB. Of course gas kinematics may show much more 
twisted features than t he stellar data (e.g. Moiseev et al .1120041: 
lEmsellem et al.ll2006|) . However, the gas in the nuclear region 
is more prone to non-gravitational forces like shocks, AGN 
jets and outflows, so may not directly probe the underlying 
gravitational potential. 

It is also worth noting that the a symmetry axis does 
not align with the secondary or the primary bar (also true 
for the sin gle bar run NB), which is consistent with what 
iMoiseev e t al. 2004 found. We do not find a clear signature in 
the hs map associated with secondary and primary bar. The 
ring in the h4 map of Figure He) is not always present for 
different projections, so it cannot be used to detect secondary 
bars. 

4.2. Slit Kinematics 

Figure [8] shows various slit profiles along the major/minor 
axis of the primary/secondary bar of run D at ^ = 405, and 
for run S. We notice that there is no central velocity disper- 
sion (a) drop in our simulated double-barred systems. So 
the cr-drop (as found in lEmsellem et al.l 1200 IL 12006) is not 
a requirement, and is not always associated with the forma- 
tion of a double-barred system. More likely a-drops are just 
the signature of newly-formed (therefore dynamically "cool") 
stars (Emsellem et al. 2001, 2006). A/^-body simulations have 
shown tha t a a-drop c an be produced in single b arred galaxies 
(fWozniak et al.ll2003l: IWozniak & Champaverti t2006), so it is 
not necessarily a unique feature of double-barred systems. 

In Figure [St a) there is a second reversal of on the major 
axis of the primary bar in the nuclear region (/? < 0.2). Such 
a hs feature is never found for a single bar case (run S) at any 
orientation. It manifests the complex asymmetric LOS veloc- 
ity distribution in the nuclear region, and is an indication of 



the decoupled secondary bar. However, the reverse argument 
is invalid: we found that such hs features are not visible at all 
orientations. 

Although our secondary bar is clearly decoupled from the 
primary bar, the maximum of the rotational velocity does not 
occur in the nuclear region (Figure [5]), instead the rotational 
velocity just rises smoothly past the region of the secondary 
bar ( the half-length of the second ary bar is a round 0.4, see 
^3.21) . This is different from Ems ellem et al.l ([2001), possi- 
bly indicating the location of the maximum of the rotational 
velocity is not a crucial factor for maintaining a nuclear bar. 

5. PATTERN SPEED DETERMINATION 

The dynamical state and evolution of barred galaxies is de- 
termined by the pattern speed of their bars. Knowledge of the 
pattern speeds of secondary bars may constrain mechanisms 
of their formation and evolution. Not much is yet known with 
certainty about secondary bar pattern speeds. Observation- 
ally, the only direct k inematic constraint onfls, ba sed on the 
Tremaine- Weinberg (Tremaine & WeinbergI 11984') method, 
was obtained for NGC 2950 by CDA03, who showed that the 
primary and nuclear bars cannot be rotating at the same rate. 
They further suggested that the nuclear bar of NGC 2950 is 
either rotating f aster or counter-r otating with respect to the 
primary bar. Maciejewski' (120061, hereafter M06), more em- 
phatically, argued that NGC 2950 has to be counter-rotating 
with respect to the primary. This would raise the prospect that 
either NGC 2950 is atypical or that counter-rotating (or possi- 
bly librating) double bars are common. However, this conclu- 
sion is based on the assumption that the Tremaine- Weinberg 
(TW) method continues to hold for nested bars, which CDA03 
suggested may not be the case. 

Our simulations provide an ideal testbed for assessing the 
reliability of measurements of fig- Here we test whether the 
simple version of the TW method as used by CDA03 using 
3 slits can recover accurately and check whether the sig- 
nature of apparent counter-rotation can occur without actual 
counter-rotation. 

5.1. The Tremaine -Weinberg method 

The TW method requires that the continuity equation is sat- 
isfied for some kinematic tracer and that the tracer's density 
can be written as T>(r^(l)-ftt). For slits parallel to the ma- 
jor axis of the disk, if V is the luminosity- weighted mean 
velocity along any such slit, and X the luminosity-weighted 
mean position along the same slit, then plotting V versus X 
results in a straight line with a slope of l^sin/. This TW 
method has been used to measure pattern spe eds in large- 
scale bars (Merrifield & Kuiiken 1995; Gersse n^t al.l Il999l: 
"Debattista et al. 2002; Asuerri et al. 2003j; iGerssen et al.l 
2003; Debattista & WilHams 2004). CDA03 showed that slits 
passing through the secondary bar did not lie on the same line 
as those passing through only the primary bar, proving that 

CDA03 argued that while can be measured from the re- 
gion outside the secondary bar, ftg cannot be obtained as eas- 
ily. Their reasons for this were two-fold: (1) disentangling the 
contribution to X and V from the primary and secondary bars 
is non -trivial and (2) the secondary bar cannot be in rigid ro- 
tation ('Louis & Ger hardlll988l:lDebattista & S"tienl l2007). vio- 
lating the assumption of the TW method that the density can 
be expressed as T^(r^(j)-ftt). Disentangling the different con- 
tributions to integrals may be possible: CDA03 presented two 
models for doing this and M06 presented another. M06 also 
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Fig. 7. — Photometrical and kinematic maps of run D and run NB. For each row from left to right are the projected surface density, mean velocities ("spider 
diagrams"), velocity dispersion, /z3, and h/[ maps, (a): Row 1, run D at r = 405 face-on. (b): Row 2, run D at r = 405 inclined at 45° with the LON of 45° relative 
to the primary bar major axis; (c): Row 3, close-up views of run D at f = 405 projected the same way as (b); (d): Row 4, as in (c) but include disk particles only; 
(e): Row 5, as in (c) but include bulge particles only; (f): Row 6, run NB at t = 400 inclined at 45° with the LON of 45° relative to the (single) bar major axis. 
The corresponding fields for run S are similar to run NB so not shown for brevity; The short and long straight line segments labels the direction of the secondary 
and primary bar, respectively (note the length of the line segment does not represent the bar length). For the projected plot, one of the dashed lines represents PA 
of the line of nodes (45°), while the other dashed is the anti-PA (135°). In the mean velocity map, the line with the connected dots shows the rough position of 
kinematic major axis, while the heavy solid curve is the zero velocity curve (kinematic minor axis). The and analyses are not preformed to bins with less 
than 100 particles. 
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t= 405.0 t= 405.0 t= 400.0 




-2-10 1 2-2-10 1 2-2-10 1 2 







Fig. 8. — Various slit profiles of run D and S. From top to bottom: surface density, mean LOS velocity, velocity dispersion, /z3, and h^. (a) left column: slit 
profiles parallel (thick curve) and perpendicular (thin curve) to the major axis of the primary bar in run D; (b) middle column: as in (a), but slits are parallel and 
perpendicular to the major axis of the secondary bar in run D. (c) right column: as in (a), but slits are parallel and perpendicular to the major axis of the (single) 
bar in run S at r = 400. 

While the disk particles support both the secondary and pri- 
mary bars, the bulge particles almost exclusively support only 
the secondary bar. Thus if we consider only bulge particles 
we have a quite clean tracer population for the nuclear bar. It 
is very unlikely that any scheme that can be devised for obser- 
vational data will ever be able to separate the secondary bar 
from the primary as cleanly as we can in our simulations. 



estimated the effect of non-rigid rotation to be less than 15%; 
this estimate was however based on the simplifying assump- 
tion that the system is 2-D. We therefore explore the effect 
of non-rigid rotation on measurements of Q.s directly with the 
simulations. The novelty of our approach lies in our ability 
to cleanly disentangle the primary and secondary bars: in our 
simulations we distinguish between disk and bulge particles. 
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5.2. Simulated TW measurements 

We therefore apply the TW-method to bulge particles only. 
We used 1 1 slits covering the full region -Ymax < Ymax, 
with Ymax = 0.3 where the nuclear bar is strongest. This cor- 
responds to slit- widths SY = 0.055 or aBi/^Y = 1 3 (in com- 
parison, the observations of CDA03 had asi/^Y = 6.4). We 
adopted an inclination / = 45° and varied the nuclear bar PA 
relative to the inclination axis, 

nuc? in the range 0° ^ '^nuc 
< 90° . We measured X and V for each slit as in .Debattistal 
(|2003l hereafter DOS): 



^ Nsiit 



iGslit 



(1) 



where Vz^i and Xi are the line-of-sight velocity and X coordi- 
nate of particle / and Nsw is the number of bulge particles in 
the slit. The sums in these definitions are over bulge particles 
only. 

To measure Vts we fit a straight line to V as a function 
of X using least- squares. As in DOS, we estimate errors on 
the slit integrals, cfx and ay, by their radial variation outside 
|X| = 0.4. We also experimented with a number of other error 
estimates including equal errors, the difference between posi- 
tive and negative Y and errors proportional to Nsw. We found 
that the most accurate measurements were obtained assuming 
weights cTy"^, which we adopt throughout. In observations 
the main uncertainties are in V and linear regression is domi- 
nated by cry, as here. We denote the slope of the fitted line as 
VLtw sin / in order to distinguish VLtw from the pattern speed, 
Vts, measured through the time evolution of the simulation. 
We quantify the typical errors in VLtw as 



An 


)-{ 











(2) 



where () represents an average over the range S0° < iprmc < 
60°, which are favorable orientations because they give large 
values of X. 

5.3. Precision of TW measurements for nuclear bars 

We start by considering the precision with which can 
be measured in the absence of a primary bar by considering 
run D at ^ = 20, before the primary bar forms but after the 
nuclear bar has saturated. Figure [9] presents the surface den- 
sity of the system; only a nuclear bar is present which is well 
traced by the bulge particles. The right panel shows the pro- 
jected surface density at / = 45° and ?/^nuc = 45° with the slits 
used superposed. The value of Vts measured from the time 
evolution is listed in Table [T] In Figure [TOl we present the TW 
measurement for the same orientation. The measured VLtw is 
accurate to better than 10%, which is the typical uncertainty 
for single bars (DOS, lO'Neill & Dubinskill200Sl) . The inte- 
grals X and V are both well-behaved and each pair of slits at 
±7 are consistent with a single straight line that matches the 
pattern speed very well. Figure fTTI summarizes the reliability 
of TW measurements for a single nuclear bar, which shows 
that Q.S can be measured to better than 10% for all reasonable 
orientations. We also experimented with a three slit configu- 
ration consisting of the central F = slit and Y = ±Foffset for 
each ^offset and found ctq increases but is still < 14%. Thus 
the TW method is well-behaved for an isolated nuclear bar. 

In Figure fT2l we show the bulge of run Dat^ = S98-415; 
the nuclear bar is prominent, and has insignificant or no elon- 
gation along the primary bar. Thus using bulge particles only 




Fig. 9. — Left: The face-on surface density (disk+bulge) in run D att = 20. 
Right: The bulge viewed at / = 45° and V^nuc =45° with the slits used in the 
TW measurement indicated by the dashed lines. 
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Fig. 10. — Left: The TW integrals as a function of slit offset Y for run D 
at f = 20 with / = 45° and V^nuc =45°. Right: the measurement of Vlj^ . The 
solid line shows O2 while the dashed line shows Vljy^ . 




'^nuc 

Fig. 11. — The fractional error in Vtjy^ at f = 20 in run D as a function of 
V^nuc- The average absolute error over 30° < V^nuc < 60° is < 10%. 

for TW measurements will result in only the nuclear bar being 
included. 

In Figure[T3]we present a single TW measurement at / = 45° 
and T/^nuc = 45° for each of the 4 times of Figure[T2l While the 
integrals themselves appear generally well-behaved, the scat- 
ter of the points about the fitted lines is larger than at ^ = 20. 
This scatter leads to a larger than in the single bar case 
for all orientations of the secondary bar, as shown in Figure 
[141 and summarized in Table [T] Other than being large at the 
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-1 -0.5 0.5 1 -1 -0.5 0.5 1 

Fig. 12. — The face-on surface density of only bulge particles in run D. 
The different snapshots are at t = 398 (top left), t = 405 (top right), t = 412 
(bottom left) and t = 415 (bottom right). The line indicates the orientation of 
the primary bar at the given time. 

smallest values of lA'l, ctq does not correlate with \X\ or |V|. 
When we fit lines to three slits as before, we find that the 
quality of the fits varies considerably (Figure [15]). We con- 
clude that observationally it is difficult to determine the uncer- 
tainty in any measurement of fts based on slit data obtained 
by CDA03 for NGC 2950. 



t 


Acj) 




crn 


20 




0.89 


0.10±0.02 


398 


135° 


0.42 


0.19±0.10 


405 


90° 


0.32 


0.10 ±0.06 


412 


45° 


0.41 


0.17±0.12 


415 


0° 


0.52 


0.09 ±0.08 



TABLE 1 

The results of TW measurements. The column Acf) gives the 

APPROXIMATE ANGLE BETWEEN THE TWO BARS. 



5.4. Interpretation 

We have demonstrated that the standard TW method on the 
secondary bar, while not wholly unreliable, is unable to re- 
cover Vis without significant uncertainty. Observationally this 
situation would be exacerbated by the need to subtract the 
contribution of the primary bar from the measured integrals, 
which we have not addressed (but see Meidt et al.l l2007). 

The amplitude of the m = 2 perturbation in bulge particles 
varies by some ±20% about the mean amplitude at all radii. 
Is the failure of the TW method for secondary bars consistent 
with the idea that non-rigid rotation leads to large errors? Ev- 
idence that this is indeed the case can be found in Figure fT4l 
which shows that the la rgest errors occur for Ac/ ) = 45° and 
= 135°. Figure 2 of IDebattista & Shed (l2007h shows that 
the amplitude of the secondary bar, A2 ~ -cos{2/S.(j)). Thus 
dA2/dt peaks at A^ = ±45°, which is in excellent agreement 
with the phases where we find the largest errors. Moreover, 
the redistribution of material being radial along the secondary 
bar, we expect that the largest errors will occur when the radial 



motions contribute more to the line-of- sight velocity. While 
some of the error at all times in Figure [14] is clearly due to 
noise, a significant part is also physical. Most importantly, 
we find that, for A^ = ±45°, the larger ijjnuc is, the larger is 
the error in VLtw- This leads us to conclude that, as argued 
by CDA03, the perturbations to the TW method due to non- 
rigid rotation are sufficiently large as to render simple mea- 
surements of Q.2 noisy at best. 

We have focused here on using slits to compare with the ob- 
servations of CDA03 and us ed only bulge particles to isolate 
the nuclear bar. Meidt et al.l (12007) present an analysis using 
an extension of the TW method which is able to disentan- 
gle multiple pattern speeds provided full 2-D velocity fields. 
They find, as here, that the pattern speed of the secondary is 
prone to larger uncertainties. However, regularization with 
that method leads to more accurate measurements of Vts- 

5.5. Comparison with NGC 2950 

Although we have shown that the TW method as used by 
CDA03 is not very accurate for secondary bars, it is not so 
grossly unreliable that we cannot consider the question of 
whether the secondary bar in NGC 2950 is counter-rotating. 
If it were, this would suggest a formation scenario for double- 
barred galaxies different from the one presented here. Sim- 
ulations have found that counter-rotating nuclear bars are 
possible if counter-rotating materia l is p resent in the disk 
(ISellwood & MerrittI Il994l: iFriedlil 1 19961: iDavies & Hunte3 
1997) and this remains a viable model if such material is 
present in a sufficiently large fraction of galaxies. 

Starting from the assumption that both bars satisfy the con- 
tinuity equation and are in rigid rotation (i.e. 11/(7? -1^/0 for 
i = s,p), CDA03 showed that the TW method for the two bars 
combined becomes 

(XpQp + XsQs) sin i = Vp-hVs = V. (3) 

The observed quantities are V and X = Xp-\- whereas the 
required quantities for determining the Qi's are A'/ = J XS/ dX 
and Vi = J Vios^i dX. Since slits can be selected to pass 
through the primary but not the secondary bar, it is possi- 
ble to derive ^}p assuming that the oscillations in the primary 
are small (in good agreement with our simulations). Ignor- 
ing the effect of non-rigid rotation, CDA03 considered two 
assumptions for A'^ in Eqn. [3] to solve for Qs in NGC 2950. 
This gave a range of possible values of 1]^, including a sec- 
ondary bar counter-rotating relative to the primary bar. Using 
the same data, M06 made a different attempt at isolating the 
secondary bar. Based on his analysis, M06 also argued that 
the secondary bar in NGC 2950 is counter-rotating. Since 
the analyses of both CDA03 and M06 ignored the non-rigid 
rotation, neither of the estimates for Qs is likely to be very 
accurate as we showed above. 

Nevertheless, we do not find in our simulations cases where 
the behavior of the integrals resembles that in NGC 2950. As 
emphasized by M06, the main characteristic of the V(Y) pro- 
file in NGC 2950 is that it becomes steeper without changing 
sign in the nuclear bar region (see Figure 3 of CDA03). This 
happens despite the fact that the two bars are on opposite sides 
of the minor axis (see Figure [J]), causing to have the oppo- 
site sign of Xp and leading to | X\ declining more rapidly in the 
secondary bar region. But instead of V also being shallower in 
this region, CDA03 found that V(Y) steepens there. For TW 
measurements of the system in Figure [2l a TW measurement 
using slits passing through both bars (now with both disk and 
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Fig. 13.— TW measurements with / = 45° and V^nuc = 45° at t = 398 (top-left), t = 405 (top-right), t = 412 (bottom-left) and t = 415 (bottom-right). The solid 
lines show while the dashed lines show Qtw- 
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Fig. 14. — The full TW analysis for only bulge particles in run D. The 
different panels are at t = 398, t = 405, t = 412 and t = 415, respectively (from 
left to right). 

bulge particles included) does not show a steeper V(Y) pro- 
file (see Figure [16]). We conclude that NGC 2950 may indeed 
have a counter-rotating primary and secondary bars. Another 
possibility might be that the nuclear bar librates about the pri- 
mary bar, which deserves more investigation in future studies. 

6. CONCLUSIONS 




Fig. 15. — The precision of TW measurements of nuclear bars in 9 models 
using 3 slits as described in the text. The different points are offset hori- 
zontally for clarity. The colors indicate t = 20: black, V^nuc — 135°: blue 
(t = 398 and t = 593), V^nuc ^ 90°: green (t = 405 and t = 599), V^nuc ^ 45°: 
yellow (t = 412 and t = 605), and V^nuc ^ 0°: red (t = 415 and t = 609). 
The error bars on the individual points show the la variations in the interval 
30° < V^nuc <60°. 



We have analyzed the photometrical and kinematical prop- 
erties of our high resolution models, and contrasted them 
when with or without a secondary bar. This study also com- 
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Fig. 16. — A TW measurement for the full system of bulge+disk att = 405 
shown in Figure [21 The solid line is for flp while the dotted line is for -fls- 
The best fit straight line to the points is shown by the dashed line. The relative 
weights of bulge particles (their "mass-to-light" ratio) have been adjusted to 
roughly reproduce A'(F) in NGC 2950. This figure is to be compared with 
Figure3ofCDA03. 

pared the simulated secondary bars with observations. 

In general the shape of secondary bars in our models is rea- 
sonable compared to observed ones. The length ratio of two 
bars, determined by various m ethods, is in the range of 0.12 t o 
0.19, in good agreement with lErwin & Sparkd (I2002L l2003h . 
We also found the overall edge-on shape of boxy bulges is 
largely unaffected by the existence of a secondary bar. At 
lower inclinations, the central density profile along the pri- 
mary bar major axis is lower than that along the minor axis, 
due to the secondary bar orienting to a different direction. 

The primary extends roughly to its corotation radius, and 
therefore fits the definition of a fast bar (see for example 
lAguerri et al.ll2 003). Although the secondary bar rotates more 
rapidly than the primary, its semi-major axis is much shorter 
than its corotation radius, even if we take the oscillation of 
the bar patterns speeds into acc ount. We did not find evi- 
dence of CR-I LR coupling (e.g. iPfenniger & NormanI 119901: 
iFriedU & Martinet 1993) in our models. 



We find that the central twist of kinematic axes is quite 
weak even if a secondary bar is present, due to the relatively 
large velocity dispersion of stars in the central region. This is 
consistent w ith the 2-D s t ellar kinematics of secondary bars 
studied in Moiseev et al." ('2004). There are no clear h4 sig- 
natures associated with the presence of secondary bars. A hs 
reversal feature may appear in the nuclear region at some fa- 
vorable orientations. We do not find a a (velocity dispersion) 
drop for our secondary bar model. It is more likely that a- 
drops are just the signature of newly-formed stars, and it is 
not necessarily a unique feature of double-barred systems. 

We showed that the Tremaine- Weinberg method is not very 
reliable even when the primary bar contribution is fully ex- 
cluded. The way in which the measurement fails is consistent 
with the proposal of CDA03, namely that the non-rigid rota- 
tion leads to internal motions that violate the stationary frame 
assumption of the method. Nonetheless, we find no exam- 
ple in our simulations where the behavior of the TW integrals 
mimics that observed in NGC 2950. Thus this galaxy may 
indeed have counter-rotating secondary and primary bars. 

The general agreement between our simulations and obser- 
vations of double barred galaxies gives us confidence that 
the simulations are capturing the same dynamics as in na- 
ture. This is especially remarkable because secondary bars 
are not merely scaled down versions of primary bars, but have 
distinctly different kinematic properties. In the absence of 
self-consistent simulations, earlier orbit-based models could 
not directly confront the challenge from observations which 
found such differences. This demonstrates the advantage of 
finally being able to simulate stellar double-barred galaxies, 
which had been puzzling for so long. 
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